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Abstract 

This work investigates the influence of low temperature and broadened Landau levels on the thermodynamic properties of two- 
dimensional electron systems. The interplay between these two physical parameters on the magnetic fleld dependence of the chemical 
potential, the specific heat and the magnetization is calculated. In the absence of a complete theory that explains the Landau level 
broadening, experimental and theoretical studies in literature perform different model calculations of this parameter. Here it is 
presented that different broadening parameters of Gaussian-shaped Landau levels cause width variations in their contributions to 
interlevel and intralevel excitations. Below a characteristic temperature, the interlevel excitations become negligible. Likewise, at 
this temperature range, the effect of the Landau level broadening vanishes. 

^Key words: two-dimensional electron systems. Landau levels, heat capacity, magnetization 
PACS: 71.10.Ca, 71.70.Da, 65.40.Ba 
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1. Introduction 

Tviro-dimensional electron systems (2DES) are widely in- 
vestigated due to their intriguing non-bulk like properties 
and their possible applications [l| . One example of interest 
is the oscillating behavior of their thermodynamic proper- 
ties with respect to a strong perpendicular magnetic field 
B at low temperature T p,l3|,|4|,|5|. This general trend is at- 
tributed to the presence of Landau levels - the discretized 
energy spectrum of a 2DES subject to a strong magnetic 
field applied normal to the 2D plane [l| . Landau levels are 
represented in the density of energy states (DOS) of a non- 
interacting 2DES as a series of Dirac-delta functions with 
gaps that are empty of electronic states [6| . 

On the side of theoretical investigations, earlier stud- 
ies have obtained Landau level broadening, in the extreme 
quantum limit, when self-consistent treatment of electron- 
impurity scattering is taken into account [7|, ^ . Another 
work, employing perturbative treatment of disorder and 
Coulomb interaction, obtained a DOS with sharp edges in 
the limit where Landau level mixing is ignored [9J. Nev- 
ertheless, reports on the DOS extracted from experiments 
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indicate broadened Landau levels. This is true whether the 
DOS is inferred from magnetization l2,| 3.llClll . heat capacity 
[4, 5] or capacitance measurements julll^]- Moreover, de- 
pending on the applied magnetic field strength, overlaps of 
Landau levels may occur which indicate existence of elec- 
tronic states in between the ideal DOS peaks. These over- 
laps have been analytically confirmed in the presence of a 
weak disorder [13|. But this derivation has no localization 
effects and it contains parameters related to the impurity 
configuration, namely, its density and distance from the 
2DES plane. Such conditions make quantitative compari- 
son with experimental results difficult. 

Others resorted to phcnomcnological calculations em- 
ploying either of these DOS shapes: Gaussian [10|, 
Lorentzian (SL 14L exponential [14| or Gaussian with a 
background [3, Isj, |5| . These depend on a single broadening 
parameter F, which according to the short range model is 
inversely proportional to B^^^ [l|. Several researchers have 
made data fits to obtain the DOS using this B dependence 
on L [2,1,11]. However, the DOS obtained using L ex B^/^ 
is found to be four times smaller than measurements in 
2DES in GaAs/AlGaAs heterostructures [13] ■ Further- 
more, the Landau level widths are predicted to oscillate 
due to the fluctuations in the screening of the disorder by 
the electrons in the 2DES la [la]- Experimental data from 
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inultiple-quantum-well structures are consistent with a F 
oscillating as a function of \/B 4]. Although there is yet 
no consensus on the i3-dependence of F, it is still advanta- 
geous to use a broadening parameter as it contains all the 
collective temperature- and magnetic field-dependence of 
the 2DES. This then makes comparison with experimental 
data of macroscopic properties easier. 

Despite the lack of a complete theory for the Landau 
level broadening caused by impurity disorder, information 
on the effects of the broadening on the thermodynamic 
quantities will provide further insights on the scattering 
mechanisms present. In this work, the behavior of the spe- 
cific heat and the magnetization as functions of F and T 
is obtained through a phenomenological approach. Here a 
Gaussian-shaped DOS is employed. This form is utilized 
since it is a commonly used model for experimental data 
and has been found to fit better than other semi-elliptical 
functions and long-ranged models [j] ■ Using various broad- 
ening parameters, the above mentioned physical quantities 
are calculated for a single sub-band in temperatures in the 
order of 0.3 K to 4.2 K and magnetic fields in the order of 
0.5 T to 15 T. 

2. The Simulation Model 

The specific heat capacity Cv and the magnetization M 
of two-dimensional electron systems are both functions of 
the magnetic field B and temperature T. The specific heat 
at constant volume of an electron gas is given as 
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where f{E, jjl, T) = l/(l-hexp [{E - fi)/kBT]) is the Fermi- 
Dirac distribution function, /i — ii{B,T) is the chemical 
potential, feg is the Boltzmann's constant, and D{B,E) is 
the density of states. 

Similarly, the magnetization for an electron gas having 
no spin splitting can be obtained from the free energy F, 



M{B, T)= - 
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where N is the electron concentration. The behavior of 
these magneto-thermodynamic properties can be deter- 
mined once the chemical potential and the density of states 
are known. 

The chemical potential /x can be obtained via a root- 
finding method for a fixed value of N . This condition is 
not a problem since experiments on 2DES usually keep A^ 
constant. For a given T and B, the value of jjl is sought 



that provides the least percentage error of not more than 
0.001% from the set value of the electron concentration in 



N = 



J fiE,^i,T)DiB,E)dE 



(4) 



In this work, N is set equal to 3.6 x 10 cm^. This value 
is chosen to be of the same order of magnitude of high 
mobility 2DES samples as found in experimental reports 
[3, 4, HE [13. 

A number of experiments support a Gaussian density of 
states D(B,E) model for two-dimensional electron systems 
[2, 3, y, |5|. The density of states at an energy E in its 
Gaussian form is given as 



D{B,E)^—Y,^^^ exp 



27rF 



{E - En) 



2r2 



(5) 



where e is the electron charge, h — hjl-K and h is Planck's 
constant. The nth Landau level is given as En = {n+-^)hujc 
where cuc = eB/m* is the cyclotron frequency. The effective 
mass TO* used here is that of carriers in a GaAs/AlGaAs 
heterostructure which is equal to 0.0667 me, where me is 
the mass of the electron. Some experimental data fits use 
Eq. ^ with an additional constant background obtained 
from a fraction of the zero- field DOS [3, |3|. This fraction 
provides an adjustable parameter dependent on experimen- 
tal factors. However, depending on B, Eq. ^ is found to 
produce Gaussian peaks with overlap between Landau lev- 
els. This overlap acts as a sort of a uniform background. 
Therefore, there is no need to introduce an additional term 
toEq. dSl). 

In our simulations, we investigate the case when the 
broadening parameter F is set to (i) fixed values, (ii) a func- 
tion of 7i?^/2 where 7 = 0.01 meV/T^/^, and (iii) an oscil- 
lating function of the filling factor 1/ similar to Ref. j^] such 
that 



T{B) ^ 0.7045^/2 + 0.296B^^^[1.8cos^{2niy) - 1] 



(6) 



where v = hN/eB. For constant F, only the results for F = 
0.2 meV and 1.0 meV will be shown. 

3. Results and Discussion 

In the absence of Landau level broadening, the chemical 
potential exhibits oscillations as a function of B with sharp 
peaks occuring at even filling factors 1^. Since the last occu- 
pied Landau level is fulUy filled for even ly, this implies that 
one needs to go to the next higher level to add one more 
electron. On the contrary, an odd v corresponds to a last 
occupied level filled to half its degeneracy. Adding an elec- 
tron will not change /x. Therefore, at odd v the chemical 
potential is a constant. 

The sharp oscillations of fJ,{B, T) are reduced and soft- 
ened upon the introduction of broadening as shown in Fig.[T] 
for T = 4.2 K at F = 0.2 nieV. A similar effect was ob- 
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Fig. 1. The chemical potential with broadened levels using F = 0.2 
meV at different temperatures. The Landau levels are located on the 
diagonal slopes of the saw-tooth like oscillations. 
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Fig. 2. The specific heat oscillations at constant T 
different broadening parameters. 
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served in 2DES in quantum wells for multiple sub-band oc- 
cupation [14]. Note, however, that despite the presence of 
broadening, [i, approaches its ideal behavior of having sharp 
peaks at T = 0.3 K. 

The Landau level broadening is reflected in the specific 
heat in Fig. [5] as spikes at the low field region. At ;/ > 4, 
the narrow peaks occur at magnetic fields when the chem- 
ical potential crosses from one Landau level to the next. 
These crossings are known as interlevel excitations [9|, ll7| . 
As B increases, the Landau levels are farther apart and it 
becomes more difficult to have interlevel excitations. This 
is reflected in the diminishing peak heights of Cy at z/ = 8 
and :^ = 6 in Fig.[2l Until finally, Cv vanishes a.\.v — 2. This 
is comparable with the results in Ref. f5| where interlevel 
excitations appear as small structures near the minima in 
the Cy oscillations at high B. 

At high i? (i/ < 4), Landau levels are far apart and the 
degree of degeneracy within a broadened level is high. Ex- 
citations brought about by an increase in B are limited 
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Fig. 3. The heat capacity for F = 0.2 meV at two different temper- 
atures. 



within the last occupied level. In this B region the intralevel 
contributions dominate ^9. Yl\. For example, Cy exhibits a 
wide peak at 2 < j/ < 4. 

This general trend of competing dominance between in- 
terlevel and intralevel contributions of Cy with B presented 
in Fig. [2] is conserved regardless of the broadening param- 
eter used. An increase in magnitude of the constant F ap- 
pears in Cy as widening of the interlevel contributions and 
narrowing of the intralevel excitations. The latter is easily 
observed in the region 2 < t^ < 4 wherein the change of Cy 
has a slower slope with B. This is consistent with the fact 
that a wider Landau level corresponds to a slower variation 
of the DOS with energy. 

When F oscillates with B one obtains a similar trend for 
Cy except for slight variations in the region where intralevel 
contributions dominate. Although not shown here, the Cy 
result for F = ^B^l"^ coincides with the result for F = 0.2 
meV at T = 4.2 K. 

The influence of the temperature on the specific heat as 
T ^ is illustrated in Fig. [3] For clarity, we will only show 
here the result for T = 4.2 K and T = 0.3 K. Reducing T 
from 4.2 K minimizes the interlevel contributions and, con- 
sequently, the magnitude of the sharp peaks at B < 4 T di- 
minishes. Until at T = 0.3 K, the effect of F vanishes. This 
shows that despite the presence of broadened Landau levels 
excitations between them are no longer possible. Only in- 
tralevel excitations contribute to Cy which is a character- 
istic of an ideal 2DES. At even v^ where i/ > 6 for F = 0.2 
meV, Cy changes from a maximum when T = 4.2K to a 
minimum when T = 0.3 K. This crossover is likewise found 
here for F = 1.0 meV at i^ > 4 and in Ref. [9| for the case 
when F ex 5^/2. 

The corresponding influence of T on M is shown in Fig.|4l 
Orbital magnetization in 2DES with electron-impurity 
scattering do not have sharp spikes with B according to 
a self-consistent theory [15{ . Moreover, experimental evi- 
dence of the de Haas- van Alphen effect displays rounding 
at the extrema of the A/ oscillations which are attributed 
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Fig. 4. The magnetization for T = 0.2 meV at T = 4.2 K for the 
dashed curve and at T = 0.3 K for the sohd curve. 
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Fig. 5. The temperature dependence of the specific heat at different 
filling factors for F = 0.2 meV. The inset shows the three Cv curves 
scaled with respect to the filling factor. 

to finite temperature and to having no states in the energy 
gaps 2^. The numerical calculations presented herein, like- 
wise, observe these rounding of M oscillations at T = 4.2 
K. Similar to the experimental results [2], the amplitude of 
the M oscillations approaches the ideal saw-tooth behavior 
with respect to _B at T = 0.3 K. The important result here 
is that these sharp oscillations are observed regardless of 
the form of the broadening, even as large as F = 1.0 meV 
where there are significant amount of states in between 
Landau levels. 

The temperature behavior of Cv corresponding to differ- 
ent fractional filling factors is illustrated in Fig. \5[ A peak 
is observed at Tp = 0.89 K. This characteristic tempera- 
ture signifies the shift from an ideal 2DES behavior to the 
regime where interlevel excitations and the Landau level 
broadening are felt. 

Below Tp, Cv has a linear T-dependence reminiscent of 
the bulk electron gas behavior [18|. Right above Tp, Cv 
decays since the intralevel contributions decrease with T 
for a given B. A similar peak was observed in Ref. [15j |. But 



they did not observe a uniform peak since their Landau 
level broadening has a non-monotonic dependence with B 
and N. Consequently, their DOS may overlap below the 
Fermi energy and have large excitation gaps beyond. The 
work presented here does not include spin splitting. It is 
worthwhile to mention that, when spin textures are present, 
Cv peak shifts to higher T with the Lande g factor for a 
given fractional v [19]. 

The same characteristic temperature is observed regard- 
less of the filling factor and the location of the chemical 
potential in the DOS. Interestingly, the Cv curves for the 
fractional v coincide if one scales the specific heat with re- 
spect to the filling factor as shown in the inset of Fig. [5l 
This Cv behavior having a peak at Tp is also observed at 
integer v. However, they do not follow this scaling behavior 
since at odd v, the chemical potential is a constant. At even 
v, the Cv peak varies in magnitude depending on whether 
the interlevel contributions dominate or not. 

4. Conclusions 

The effect of Landau level broadening and temperature 
is considered in two-dimensional electron systems using a 
Gaussian-shaped density of states. The broadening param- 
eters considered took different forms such as different con- 
stant values, a square-root dependence on the magnetic 
field and an oscillating function with respect to the filling 
factor. The observed B behavior of the chemical potential, 
specific heat and magnetization show the same trend re- 
gardless of the form of F. The broadening only plays a role 
in the width of the oscillations with B for each interlevel 
and intralevel contribution. A characteristic temperature 
is observed such that the 2DES approaches its ideal elec- 
tron gas behavior. Below this characteristic temperature, 
the Landau level broadening no longer affects the behavior 
of the thermodynamics properties. This indicates that scat- 
tering mechanisms present in 2DES are suppressed below 
Tp. The magnitude of Tp is obtained here using the electron 
concentration and effective mass value for a GaAs system. 
It would be interesting to know if experiments are able to 
measure Tp. If not, this should be a ground for theoretical 
studies to find out why Tp is masked in experiments. 
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